Emergence of fractal behavior in condensation-driven aggregation 
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We investigate a model in which an ensemble of chemically identical Brownian particles are con- 
tinuously growing by condensation and at the same time undergo irreversible aggregation whenever 
two particles come into contact upon collision. We solved the model exactly by using scaling theory 
for the case whereby a particle, say of size x, grows by an amount ax over the time it takes to 
collide with another particle of any size. It is shown that the particle size spectra of such system 
exhibit transition to dynamic scaling c(x,i) ~ t~^ '<j>(x /t z ) accompanied by the emergence of fractal 
of dimension df = 1 + 2a ■ One of the remarkable feature of this model is that it is governed by a 
non-trivial conservation law, namely, the df moment of c(x, t) is time invariant regardless of the 
choice of the initial conditions. The reason why it remains conserved is explained by using a sim- 
ple dimensional analysis. We show that the scaling exponents /3 and z are locked with the fractal 
dimension df via a generalized scaling relation j3 = (1 + df)z. 
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PACS numbers: 61.43.Hv, 64.60.Ht, 68.03.Fg, 82.70Dd 



I. INTRODUCTION 

The formation of clusters by aggregation of particle, its 
underlying causes and consequences, is one of the most 
fundamental yet challenging problem of many processes 
in physics, chemistry, biology and engineering. Exam- 
ples include aggregation of colloidal or aerosol particles 
suspended in liquid or gas 0, 0, 0| , polymerization [3] , 
antigen-antibody aggregation 5] and cluster formation 
in galaxy Q . Such a wide variety of applications has re- 
sulted in numerous studies focusing mostly on the kinetic 
and geometric aspects of the problem. The kinetic aspect 
is well studied and well understood by theory, experiment 
and by numerical simulation. The first successful theo- 
retical model was proposed more than one hundred years 
ago by von Smoluchowski and it still remains the only 
analytical model which has provided much of our theo- 
retical understanding Q ■ The definition of this model is 
trivially simple. It is assumed that initially an ensem- 
ble of chemically identical particles undergo sequential 
aggregation upon collision. 

The kinetics of aggregation by the Smoluchowski equa- 
tion has been extensively studied in and around the 1980s 
and significant contribution was made during this period 
especially on the scaling theory and gelation transition 
[1, [^, [To| . One of the most striking results is that the 
concentration c(x, t) of particles of size x at time t ex- 
hibits dynamic scaling 



c(x,t) ~ s(t)~ e (j>(x/s{t) 



(1) 



in the limit t — > oo, where is the scaling function, 
s(t) is the mean particle size and the conservation of mass 
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principle tunes the mass exponent to an integer value 
= 2^, [ll| . The structure of the above scaling ansatz 
is highly instructive as it has been found in many seem- 
ingly unrelated phenomena. It implies that there must 
exists a common underlying mechanism for which such 
disp arate systems behave in a remarkably similar fashion 
On the other hand, the insights into the geometric 
aspect is mostly provided by experiments and numerical 
simulations and these studies reveal that when particles 
agg regate almost always scale-invariant fractals emerge 
Unfortunately, there does not yet exist an ana- 
lytically solvable model which could help us know why 
fractals are ubiquitous in the aggregation process. 

In addition to aggregation, particles may also grow in 
size by condensation, deposition or by accretion. For 
instance, aggregation in vapor phase or in damp envi- 
ronment particles or droplets may continuously grow by 
condensation [l3|, 0, [l5|, [l(| . It is also well-known that 
aerosol or colloidal particles are often not stable rather 
their sizes may evolve via aggregation and condensation 
leading to gas-to-particle conversion. However, in the 
absence of impurity such as dirt or mist, the condensa- 
tion can only take place on the existing particles without 
forming new nuclei provided the concentration of parti- 
cles present is sufficiently high and the supersaturation 
is sufficiently low [l|, [H[ . This type of growth is known 
as the heterogeneous condensation. To this end, we re- 
cently proposed a simple condensation-driven aggrega- 
tion (CDA) model and discussed the kinetic aspect of 
the problem [_t| . In this article, we present an alterna- 
tive method to solve the CDA model and kept our focus 
mainly on its geometric aspects instead. We show ana- 
lytically that the resulting system can be best described 
as fractal and quantified by its dimension df which de- 
creases with increasing strength of growth by condensa- 
tion. Interestingly, we find that the key results of the 
CDA model are connected, in one way or another, to 
the fractal dimension df. For instance, the dfth mo- 
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ment of the distribution function c(x, t) is a conserved 

i 

quantity, the mean particle size grows with time as , 
in terms df we can write a generalized scaling relation 
/3 = (1 + df)z, etc. To test our analytical predictions, 
we have performed extensive numerical simulation and 
found that analytical results are in perfect agreement 
with numerical data. 

The rest of the paper is organized as follows. In Sec. 
II, we give a detailed description of the CDA model in- 
cluding its algorithm. In Sec. Ill, some of the key fea- 
tures of the model are discussed. In Sec. IV, we give 
a simple dimensional analysis to the governing equation 
of the CDA model in an attempt to gain deeper insight 
into the scaling theory. We applied the scaling theory in 
Sec. V to obtain the solution for the distribution function 
c(x,t). In Sec. VI, we invoke the idea of fractal analysis 
to the CDA model. Finally, in Sec. VII we discuss and 
summarize our work. 



II. THE MODEL 

Chemically identical particles in aggregation process 
are typically characterized by their mass or size and 
shape. However, if the particles are one dimensional then 
size or mass is the only dynamical variable. Therefore, 
within a given class of units of measurement both mass 
and size can be described by the same numerical value as 
they differ only by a proportionality constant. However, 
this is not true in the case of higher dimensional particles. 
The Smoluchowski model is inherently one dimensional 
and hence size and mass can be used interchangeably. In 
the CDA model, we characterize each particle by the size 
it assumes upon aggregation till it takes part in further 
aggregation. The extent of growth by condensation can 
be quantified by the growth velocity defined as the ratio 
of the net growth and the elapsed time during which this 
growth occurs. The most natural choice for the elapsed 
time is definitely the collision time. The growth velocity 
is then fully specified if we know the amount of growth of 
a given particle which occurs between collisions. For this, 
we assume that the net growth of a particle between col- 
lisions, in the most generic case, is directly proportional 
to the size by which it is characterized. That is, a parti- 
cle which is just born upon aggregation with size x will 
have its size equal to x + ax whenever it collides with 
another particle regardless of the amount of time it takes 
to collide. This is not at all a bad assumption since such 
a choice makes the growth velocity stochastic in nature 
as the growth size and the collision time both become 
random in character. 

For numerical simulation, one may think of keeping 
a logbook where the sizes of the particles are registered 
each time they take part in aggregation. Initially, sizes 
of all the particles in the system are registered in the 
logbook. The rules these particles then have to follow at 
each step during simulation are: 



1,1,1,1,1,1,1,1,1,1 



4,1,1,1,1,1,1,1,1 



10,1,1,1,1,1,1,1 



22,1,1,1,1,1,1 



10,4,1,1,1,1,1 



16,1,1,1,1,1,1 



10,4,1,1,1,1,1 



FIG. 1: Schematic representation of the model for a = 1 is 
given by using monodisperse initial condition as an example. 



(i) Two particles are picked randomly from the system 
to mimic random collision via Brownian motion. 

(ii) The sizes of the two particles are increased by a 
fraction a of their respective sizes in the logbook 
to mimic the growth by condensation. 

(iii) Their sizes are combined to form one particle to 
mimic the aggregation process. 

(iv) The logbook is updated by registering the size of 
the new particle in it and at the same time deleting 
the sizes of its constituents from it. 

(v) The steps (i)-(iv) are repeated ad infinitum to 
mimic the time evolution. 

In order to illustrate how these rules of the model work 
for monodisperse initial condition, we give a simple ex- 
ample in Fig. (I) using an evolutionary tree based ap- 
proach. 

The CDA model can also be understood by a reaction 
scheme written as 



v(x t t) 



A x {t) + A y (t) -A A (a+1){x+y) (t + T), 



(2) 



where A x (t) denotes the aggregate of size x at time t 
and t is the elapsed time. This reaction scheme can 
be described by the following generalized Smoluchowski 
(GS) equation 

r d d 1 f°° 

+ ~ v ( x ^) c(x,t) = -c(x,t) / K(x,y)c(y,t)dy 
J Jo 



dt dx 
I 



2 Jo 



dyK(y, x - y)c(y, t)c(x -y,t). 



(3) 



The second term on the left hand side of the above equa- 
tion accounts for the growth by condensation with veloc- 
ity v(x, t). On the other hand, the first (second) term on 
the right hand side of Eq. ([3]) describes the loss (gain) 
of size x due to merging of size x {{x — y)) with particle 
of size y. However, the GS equation can only describe 
the CDA model if the growth velocity v(x,t), the colli- 
sion time r, and the kernel K(x, y) are suitably chosen 
as required by the rules (i)-(v). For instance, according 
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FIG. 2: Plots of ln(s) versus ln(t) are shown for three different 
a values but w ith the same monodisperse initial conditions 
in each case. The lines have slopes equal to (1 + 2a) revealing 
the same growth law as predicted by Eq. (I10|l . 



to rule (ii) of our model, the net growth of a particle of 
size x between collisions is ax. To obtain a suitable ex- 
pression for the elapsed time we do a simple dimensional 
analysis in Eq. and immediately find that the inverse 
of J °° K(x,y)c(y,t)dy is the collision time t(x) during 
which the growth ax takes place fl8j . The mean growth 
velocity between collisions therefore is 



v{x, t) 



ax 
t(x) 



<>.x 1 dyK(x,y)c(y,t). 
'0 



(4) 



The rule (i) on the other hand says that a given parti- 
cle can collide with any particle in the system with an 
equal a priori probability regardless of their size. This 
can be ensured only if we choose an aggregation kernel 
independent of its argument and hence we set 



K(x,y) = 2, 



(5) 



for convenience. 



III. SOME OF THE BASIC FEATURES 

It is noteworthy to mention that the distribution func- 
tion c(x,t) itself is not a directly observable quantity but 
its various moments are. Therefore, one often finds it 
more convenient to deal with its moment than the func- 
tion itself. The fcth moment of c(x, t) is defined as 



Mfc(i) = / x k c(x,t)dx, with k > 0. 



(6) 



Incorporating it in Eq. §5§ after substituting Eqs. ((31) 
and ([5]) in it we obtain 



dM k (t) 
dt 



OO pOO 



dxdyc{x, £)c(y, t) 



(7) 



x (x + y) k + (ak - l)(x k + y k ) 

In the case of a — which describes the classical Smolu- 
chowski (CS) equation we find dM ^^ = and hence 



FIG. 3: Plots of ln(s) versus ln(t) for monodisperse (all the 
particles are chosen to be of unit size) and polydisperse (e.g., 
the numerical value for the size of the 500 particles are cho- 
sen randomly from the interval [1, 10, 000]) initial conditions. 
Two parallel lines prove that the growth-law for s(t) given by 
Eq. (llOjl is independent of initial conditions. 

Mi(t) = Jg 00 xc(x,t)dx is a conserved quantity . This is 
well known as the conservation of mass principle. Ob- 
viously, this principle is no longer obeyed in the CDA 
model because of the growth by heterogeneous conden- 
sation. To check we solve Eq. ([7]) for Mi(t) = C(t) and 
find that 



C(t) ~ 1 2 



(8) 



in the long time limit. Hence, the growth of the total 
mass or size exhibits non-universal behaviour in the sense 
that its exponent depends on the parameter a. It con- 
firms that the conservation of mass principle is violated 
Va > 0. Then the question remains: Can there be an- 
other conservation law which the system should obey as 
it evolves? One cannot find a straightforward answer to 
this question from Eq. ([7]) except the a — case. We 
now solve Eq. (J7J for k = to obtain the solution for the 
total number of particles Mo(t) = N(t) present at time t 
and find that in the long time limit it decays algebraically 
with a universal exponent 



N(t) ~ r 



(9) 



In other words, the number density N, evolves following 
the same differential equation as the one we would obtain 
for the CS equation. It therefore confirms that conden- 
sation takes place only on the existing particles without 
forming new nuclei. 

Using the solutions for the first two moments in the 
definition for the mean particle size s(t) = C(t)/N{t) we 
obtain the following growth-law 



s(t) ~ t 



l+2a 



(10) 



in the long time limit. To verify this, we performed nu- 
merical simulation based on the rules (i)-(v). However, 
in order to manipulate the numerical data we define time 
t = 1 /N since the number of particles present in the sys- 
tem determines how fast or slowly the aggregation pro- 
cess should proceed. We first use the monodisperse ini- 
tial condition where all the particles are assumed to be 
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of unit size. In Fig. (2), we present plots of ln(s) versus 
ln(t) from the resulting data and find three straight lines 
for three different values of a. The slopes of these lines 
satisfy the relation z = 1 + 2a, V a > which clearly 
shows algebraic growth of s(t) as predicted by Eq. (fTO)) . 
The next important thing is to check if the initial distri- 
bution of particle size has any effect in this growth law. 
To find this out we simulated the model for several differ- 
ent polydisperse initial conditions and collected data for 
the mean particle size s(t) against time t. In one of the 
instances, we picked initially 500 particles of size chosen 
randomly from the interval [1,10,000] and let the pro- 
gram run in the computer following the rules (i)-(v) of 
the algorithm. In Fig. (3), we again present plots of ln(s) 
versus ln(£) from the resulting numerical data and put it 
together with the corresponding plot for the monodis- 
perse initial condition to see the contrast. Surprisingly, 
we find two parallel lines which clearly implies that the 
exponent of the growth law is universal in the sense that 
it is independent of the initial conditions. That is, one 
may choose any number of his choice to replace the parti- 
cles characterized by 1 in the first box [t — 0) of Fig. (1) 
and let them play following the rules (i)-(v). One would 
still obtain the same growth-law for the mean particle 
size as the one for the monodisperse initial condition. 



IV. DIMENSIONAL ANALYSIS 

There are two governing parameters x and t in the GS 
equation. However, according to Eq. (fT0|) the size of the 
particle can be expressed in terms of time. Therefore, 
only one of the variable, say, time t can be taken as an 
independent parameter. The other governing parameter 
such as x and the governed parameter c(x,t) both can 
be expressed as a function of t alone. We already know 
from Eq. iflU)) that t z with kinetic exponent 



z = 1 + 2a, 

have the dimension of length [s(t)] 
define a dimensionless quantity 



(11) 

L and hence we can 



(12) 



On the other hand, applying the power-monomial law 
for the dimension of physical quantity we can write a di- 
mensional relation c(x, t) ~ where the exponent j3 
assumes such a value that makes i - * 3 bear the dimen- 
sion o£c(x,t) [H. We can therefore define yet another 
dimensionless quantity <j> as follows 



c(x, t) 



(13) 



Now, within a given class one can pass from one units of 
measurement to another system of units of measurement 
by changing t by an arbitrary factor leaving the other 
factor unchanged. Upon such a transition the quantity 



on the right hand side of Eq. (| 13|) remain unchanged 
since the left hand side is a dimensionless quantity. It 
means that the quantity c ^jJ) vis-a-vis (j) can only be a 
function of another dimensionless quantity and the only 
dimensionless governing parameter is £ given by Eq. (|12[) . 
We therefore find that Eq. ([L3")) leads to the following 
scaling ansatz 



(14) 



or can be expressed in the form of Eq. ([T]) if we set 
j3 = 9z and use Eq. (TIT)]) thereafter. Existence of scaling 
means the following. The quantity c(x, t) that depends 
on two variables x and t is considered to admit scaling 
if the two variables combine into one variable such away 
that it can be expressed as Eq. (fT4")l . The fact that two 
variables combine into one variable leads to an enormous 
simplification in finding the solution to the problem as 
we shall see below. 



V. SCALING SOLUTION 

To check if the solution of the GS equation exhibits 
scaling or not, we substitute Eq. (fT4")) together with [3 = 
9z in the GS equation after substituting Eqs. ([4]) and ([5]) 
in it and obtain 



t" 



where, 



1 



no 

<t>{v)<Kt 



(2(1 + a)$ ^(0 + 2a$ £ 



r))drj 



(15) 



F(0 = 9z4>(0 + z^'(0, 



(16) 

and $o = Jq 00 <X£)rf£ is the zeroth moment of the scaling 
function Note that the right hand side of Eq. ([15)) 

is dimensionless while the left hand side is not. Thus the 
dimensional consistency requires 



1 



- 1 



(17) 



Substituting the z value from Eq. pip in the above equa- 
tion we obtain the mass exponent 



2 + 2a 
l + 2a' 



(18) 



The values for the exponents 9 and z are exactly the same 
as obtained by the exactly solvable method, namely the 
Laplace transformation and the method of characteris- 
tics, in Ref. [13] • To obtain the complete scaling or 
self-similar solution of the GS equation we still have to 
find </>(£). For that, we substitute the value of z and 9 in 
Eq. ([]]§ to get 



1 + 2a(l - * ) + / ~ V)dv 

J d£ Jo 

+ [2(l + a)(l-$ )"U(0 = (19) 



5 



The solution of the problem thus reduces to finding the 
solution of an ordinary integro-differential equation for 
the scaling function </>(£). 

We now multiply on both sides of Eq. (TH))) by £™ and 
integrate from £ = to £ = oo to obtain an equation for 
the nth moment <E>„ of (f>(£) (or Mellin transform <I>(n + 
1) = $„) which can be written in the closed form 

n 

2on(l - $o)$„ + (n - 1)$„ + 2*„$ = ^ n Cr*r*n-r, 

r=0 

for integer value of n only. By setting n — in the above 
equation we find 



$ ($o - 1) = 0, 



(21) 



and the only non-trivial solution of this equation is $q = 
1. Using it back in Eq. I|19p gives 



(22) 



To find a solution of this equation we set $o = 1 in 
Eq. ((20)) and find the following hierarchy of equations 
for different integer n values and a few of these are 



$ 4 = I[ 8 $3$ 1+6 $2]. e t C j 

The solution $n = 1 and its definition 



(23) 



implies 



*o = / em<%, 

Jo 



m = e 



(24) 



(25) 



which is in fact the inverse Mellin transform of <E>(1) 
$o- It can also be verified by substituting it in Eq. 
Indeed, one can check that this solution does satisfy the 
hierarchy of all the relations in Eq. (|23p and simultane- 
ously it solves Eq. (f2"2"|) . 

Substituting Eqs. {TUH, (H) and © in Eq. Q we 
can finally obtain the scaling solution of the GS equation 



(26) 



This is exactly what we found in Ref. from the ex- 
plicit time dependent solution by using the limit t — > oo. 
However, the exact solution in [17j | was obtained only for 
the monodisperse initial condition. The advantage of us- 
ing the scaling theory is that we do not need to specify 
the initial condition. It implies that the various essential 
features or laws of the CDA model should remain inde- 
pendent of the initial condition provided the number of 
particles present at time t = is sufficiently large and 
one allows the system to run for sufficiently long time. 




FIG. 4: Plots of the ( i+2a ) ^ moment of the particle size dis- 
tribution function versus time for different initial conditions 
(monodisperse and polydisperse) and for different a values are 
drawn using numerical data. A set of horizontal lines clearly 
prove that the generalized conservation law is always obeyed 
regardless of the initial conditions or the a value. 



We can now obtain the solution of Eq. ([7]) , the kth mo- 
ment of c(x, t) by using Eq. (|26[) in its definition given 
by Eq. ^ to give 



M k (t) 



(27) 



A surprising feature of this solution is that it implies the 
existence of a non-trivial conservation law namely the qth 



moment is independent of time if q = 



l+2a 



In numerical 



simulation, it means that the sum of the ( — ^ j th power 

of all the particles is a constant as a function of time. To 
check this, we plot ln(M 9 ) with q = 1+ 1 2a against ln(i) 
from the numerical data for different a values and for 
different initial conditions. We find a set of parallel hor- 
izontal lines for all a > values (see Fig. (4)) regardless 
of initial conditions value revealing that our analytical 
result is in perfect agreement with the numerical simula- 
tion. One may ask: What is so special about this q value 
that makes this moment a conserved quantity? To find 
an answer to this question we invoke the idea of fractal 
analysis which is discussed below. 



VI. FRACTAL ANALYSIS OF THE CDA 
MODEL 

In this section, we intend to address the geometric as- 
pect of the CDA model by invoking the idea of fractal 
analysis. Before doing so we find it worthwhile to appre- 
ciate the following. Theoretically, particles in the CDA 
model can be considered to be inside a row of boxes form- 
ing a Id lattice with each box as one lattice point. Ini- 
tially, each of these boxes contains one particle character- 
ized by one number whose size distribution depends on 
the initial particle size distribution. In this sense, parti- 
cles are in fact embedded in a space of dimension equal to 
one. The geometry of the resulting system therefore will 
be called fractal if the dimension of the distribution of 
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particles is less than one and greater than zero. To know 
exactly what this value is we define C(t) as the measure 
which is the sum of all the aggregates in different boxes 
at time t and obviously it is an ever growing quantity 
against t according to Eq. ([8]). 

To quantify the measure C, one can use a suitable yard- 
stick and find an integer number N needed to cover the 
measure C. The most suitable candidate for the yardstick 
in the context of the CDA model is the mean particle 
size s(t) which will always give the number N an integer 
value. That is, the size of the measure C(t) can be quan- 
tified by the number N(s). Using k — in Eq. (f2"T)) . we 
can easily see that the number N(s), when expressed in 
terms of s(t), exhibits power-law 



N{s) 



with exponent 



df = 



1 



l + 2a 



(28) 



(29) 



which is highly significant for the following reason. Note 
that when the number N is obtained by measuring a 
given measure with a suitable yardstick and find that 
it exhibits power-law against the size of the yardstick, 
then the exponent of the power-law is widely known as 
the Hausdorff-Besicovitch (H-B) dimension [21|. On the 
other hand, the H-B dimension is called fractal if it is 
non-integer and at the same time if it is less than the 
dimension of the embedding space. It implies that the 
exponent df of Eq. (f2"5| is the fractal dimension of the 
resulting system since, according to Eq. (|2"9"j). df is not 
only non-integer V a > but also less than the dimen- 
sion of the embedding space d = 1. It is noteworthy that 
the size of the fractal that emerges in the CDA model is 
continuously growing with time but at the same time it 
preserves its dimension which is ensured by the conserva- 
tion law. Within the rate equation approach, such frac- 
tal analysis was first done by Ben-Nairn and Krapivsky 
in the context of the stochastic Cantor set [22j and later 
one of us applied it successfully in several different sys- 
tems (23[. To verify our analytical result, we have drawn 
ln(iV) versus ln(s) in Fig. (5) from the numerical data 
collected for a fixed initial condition but varying only the 
a value. On the other hand, in Fig. (6) we have drawn 
the same plots for a fixed a value but varying only ini- 
tial conditions (monodisperse and polydisperse) . Both 
figures show an excellent power-law fit as predicted by 
Eq. (|28p with exponent exactly equal to df regardless 
of the choice we make for the initial size distribution of 
particles in the system. 

We shall now show that the various interesting results 
of the CDA model can be expressed in terms of the fractal 
dimension df. For instance, we can use the expression for 
the fractal dimension df in Eq. (jTTj) to obtain z = 

Using it in Eq. (fT0|) we obtain the following growth law 
for the mean particle size 




x=0.25 □ 
<x=O.B O 
ct=1.0 Z\ 



FIG. 5: Plots of ln(iV) vs ln(s) are drawn for three different 
a values keeping the same initial condition. The lines have 
slopes equal to 1 ^ 2a which is exactly what was predicted by 
the theory. 





FIG. 6: Two parallel lines resulting from the plots of ln(JV) 
versus m(s) for monodisperse and polydisperse initial condi- 
tion reveal that N ~ s~ d f is independent of initial conditions. 



We find it worthwhile to mention here, as a passing note, 
that a similar growth law has also been found experimen- 
tally by Weitz et al while studying the diffusion-limited 
cluster-cluster aggregation [24j]. Also, the mass exponent 
9 can be expressed in terms of df by using Eq. (|29p in 
Eq. (JTHJl to give 



= 1 



'I- 



(31) 



We thus find that 8 always satisfies the inequality 9 < 
2, V a > and the inequality becomes equality 8 = 2 
only if a = which corresponds to CS model. One can 
interprete the above expression for the mass exponent 8 
as the sum of the fractal dimension df and that of its 
embedding space 1. Using Eq. (f3"Tj) in j3 = 9z we can 
further write a generalized scaling relation 



(32) 



s(t) ~ t d f 



(30) 



It is interesting that a similar expression for the expo- 
nents 9 and (3 have also been found in other phenomena 
which indicates that these results are universal in char- 
acter m . 

To further support our results, we once again use 
the simple dimensional analysis. According to Eq. dTJ 
the physical dimension of c(x,t) is [c] = since 
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[s(t)} = L and 6 = 1 + df. On the other hand, the con- 
centration c(x, i) is defined as the number of particles per 
unit volume of embedding space (V ~ L d where d = 1) 
per unit mass (M) and hence [c] = L _1 M~ l . Now ap- 
plying the principle of equivalence we obtain 

M(L) ~ L d/ . (33) 

This relation is often regarded as the hallmark for the 
emergence of fractality. An object whose mass-length 
relation satisfies Eq. (|33|) with non-integer exponent is 
said to be fractal in the sense that if the linear dimen- 
sion of the object is increased by a factor of L the mass 
of the object is not increased by the same factor. That 
is, the distribution of mass in the object becomes less 
dense at larger length scale. It proves that the splits of 
the mass exponent into dimension of the fractal (df) and 
that of its embedding space (d = 1) is consistent with 
the definition of the distribution function c(x, t) as well. 
It is interesting to note that such a simple dimensional 
analysis can also provide us with an answer to the ques- 
tion: Why is the moment Md f = J Q x df c(x,t)dx a con- 
served quantity? For an asnwer, we find it conventient to 
look into the physical dimension of its differential quan- 
tity dMd f = x df c(x,t)dx. Using the physical dimension 
[x] = L and [c(x, i)] = in the expression for 

dMd f , we immediately find that it bears no dimension 
and so is the quantity Md f . Recall that the numeri- 
cal value of a dimensionless quantity always remain un- 
changed upon transition from one unit of measurement to 
another within a given class. In the context of the CDA 
model it implies that the numerical value of Mdf remains 
the same despite the fact that the system size continues 
to grow with time. It is due to this reason that we find 
that the d/th moment of c(x,i) is a conserved quantity. 
We thus see that the simple dimensional analysis proved 
to be very useful in gaining the comprehensive explana- 
tions of various results of the CDA model which we have 
been longing for. 

VII. DISCUSSION AND SUMMARY 

In this work, we studied the geometric aspect of the 
condensation-driven aggregation model that we recently 
proposed. In the present work we first gave a simple 
dimensional analysis to the generalized Smoluchowski 
equation as we found it provided not only a deeper in- 



sight but also, at the same time, an elegant way to look 
into the problem. We then applied the scaling theory 
and shown that the GS equation admits simple scaling 
only if z(6 — 1) = 1. That is, the solution for parti- 
cle size spectra exhibits transition to dynamic scaling 
c(x,t) ~ t~ 6z (j>(x/t z ) with scaling function </>(£) ~ e~^. 
Substituting the solution for the distribution function 
c(x, t) into the definition of the nth moment shows that 
the moment of order equal to 1+ 1 2a is a conserved quan- 
tity throughout. In an attempt to know exactly why this 
value is so special we invoked the idea of fractal analysis 
and found that it is in fact the value of the fractal dimen- 
sion of the resulting system. The expression for fractal 
dimension df = 1 ^ 2a states that as the extent of growth 
by condensation increases the dimension of the measure 
or the object decreases which is quite counter intuitive. 

To summarize, we found that the fractal dimension df 
plays a pivotal role in describing and understanding the 
geometric aspect of the CDA model. For instance, the 
dynamics of the system is governed by a conservation 
law which is the d^ 1 moment of the distribution func- 
tion c(x,i). The exponent of the algebraic growth-law 
for the mean particle size is equal to -jj. In terms of df 
we can express the mass exponent and the scaling rela- 
tion in their generalized form such as = 1 + df and 
(3 = (1 + df)z respectively. A simple dimensional anal- 
ysis to the distribution function c(x, t) and the use of 
= 1 + df led to the well known mass-length relation 
M(L) ~ L df . The dimensional analysis also revealed 
that the mass exponent 9 is in fact equal to the sum 
of the fractal dimension df and that of its space where 
it is embedded. Besides, we have shown that the d/th 
moment Md s is actually a dimensionless quantity and 
we argued that this is exactly the reason why Md f re- 
mains time invariant. We thus see that appreciation of 
the exponent d/ as the fractal dimension has provided 
a self-consistent explanation to all the results which are 
found to be independent of the initial particle size dis- 
tribution. Moreover, we have shown that the results are 
independent of initial conditions. The ideas developed 
in this paper could be taken further by investigating the 
CDA model for aggregation kernel K(x, y) — (xy) w . This 
would be an ideal case in order to study how the onset of 
gelation is modified, if at all, by the presence of growth 
by heterogeneous condensation. We hope to address this 
issue in our future endeavour. 
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